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,__, ; Abstract 

In a recent paper Bender and Mannheim showed that the unequal-frequency fourth-order deriva- 
t— i ■ tive Pais-Uhlenbeck oscillator model has a realization in which the energy eigenvalues are real and 

^-; 

_J , bounded below, the Hilbert-space inner product is positive definite, and time evolution is unitary. 

Central to that analysis was the recognition that the Hamiltonian Hp\j of the model is VT symmet- 



o 



ric. This Hamiltonian was mapped to a conventional Dirac-Hermitian Hamiltonian via a similarity 



transformation whose form was found exactly. The present paper explores the equal-frequency limit 
of the same model. It is shown that in this limit the similarity transform that was used for the 
unequal-frequency case becomes singular and that Hp\j becomes a Jordan-block operator, which is 
nondiagonalizable and has fewer energy eigenstates than eigenvalues. Such a Hamiltonian has no 
Hermitian counterpart. Thus, the equal- frequency VT theory emerges as a distinct realization of 
quantum mechanics. The quantum mechanics associated with this Jordan-block Hamiltonian can 
be treated exactly. It is shown that the Hilbert space is complete with a set of nonstationary solu- 
tions to the Schrodinger equation replacing the missing stationary ones. These nonstationary states 
are needed to establish that the Jordan-block Hamiltonian of the equal-frequency Pais-Uhlenbeck 
model generates unitary time evolution. 



I. INTRODUCTION 

A decade ago it was discovered that the non-Dirac-Hermitian Hamiltonian H = p 2 + ix 3 
has an entirely real quantum-mechanical energy spectrum [l|. The reason for this unexpected 
spectral reality is that this Hamiltonian has an underlying VT symmetry and that this 
symmetry is unbroken; that is, the energy eigenstates are also eigenstates of the VT operator. 
In general, whenever a Hamiltonian has an unbroken VT symmetry, its energy spectrum 
is real [2|, |3j. Unbroken VT invariance serves as an alternative to Dirac Hermiticity in 
quantum theory. Moreover, for any non-Hermitian Hamiltonian that has a complete basis 
of real-energy eigenstates, there necessarily exists a similarity transformation that brings it 
to diagonal Hermitian form. Because the transformation is a similarity rather than a unitary 
one, such unbroken-'PT-symmetric Hamiltonians are really Dirac- Hermitian Hamiltonians 
written in a disguised form in terms of a skew basis. However, even if a non-Hermitian 
Hamiltonian has a real eigenspectrum, it is not automatically diagonalizable because it 
might be of Jordan-block form. In such an event no similarity transformation exists and the 
VT sector then emerges as a distinct and self-contained realization of quantum mechanics. 

In this paper we study an exactly solvable model, namely, the equal-frequency version of 
the Pais-Uhlenbeck oscillator [4| in which a distinct VT realization arises and for which the 
■pT-symmetric Hamiltonian has no Dirac-Hermitian counterpart. We explore the quantum 
mechanics associated with such a self-contained VT quantum realization and give special 
attention to the lack of completeness of the energy eigenstates that is characteristic of 
Jordan-block Hamiltonians. Previously, we were motivated to study the unequal-frequency 
version of the Pais-Uhlenbeck model because of its ghost problem [5[. Here, we study the 
Jordan-block structure of the equal-frequency model by extending the techniques that were 
developed in Ref. [5| to resolve the ghost problem in the unequal-frequency version of the 
model. We show that even though the Pais-Uhlenbeck-model Hamiltonian develops Jordan- 
block structure in the equal- frequency limit, the unitarity of the theory is not lost. 

This paper is organized as follows: Section [III provides a brief summary of VT quan- 
tum mechanics and Sec. IIHI reviews the results of Ref. [5] for the unequal-frequency Pais- 
Uhlenbeck model. In Sec. IIVI we construct the Fock space associated with the unequal- 
frequency Pais-Uhlenbeck model. Then, in Sec. [V]we use the results of Sec. [IV] to construct 
the equal-frequency Fock space, and we show that a Jordan-block form for the Hamiltonian 



with its incomplete set of energy eigenstates arises in the equal-frequency limit. In Sec. |VT] 
we construct the eigenfunctions of the unequal-frequency theory, and in Sec. IVHI we take 
their limit to construct the eigenfunctions of the equal-frequency theory. In doing so, we 
discover what happened to the missing energy eigenstates. In Sec. IVIIII we present some 
conclusions and comments. Finally, in the Appendix we discuss the implementation of the 
Lehmann representation in theories having non-Hermitian Hamiltonians. 

II. BRIEF SUMMARY OF VT QUANTUM MECHANICS 

There has been much research during the past few years on VT quantum mechanics; early 
references include Refs. |l|, |2|, [6| and some recent reviews may be found in Refs. |3j, ID]- A 
VT quantum theory is one whose dynamics is governed by a Hamiltonian H that commutes 
with the VT operator. (Here, V is the parity operator, which performs spatial reflection, 
and T is the time- reversal operator.) The appeal of VT quantum mechanics is that even if 
a Hamiltonian is not Dirac Hermitian it will still have real energy eigenvalues whenever it 
is VT symmetric and all of its eigenstates are also eigenstates of the VT operator. 

For diagonalizable "PT-symmetric Hamiltonians it is convenient to construct an operator 
C, which obeys the three simultaneous algebraic operator equations 

C 2 = l, [C,VT]=0, [C,H]=Q. (1) 

The first two of these equations are kinematical, while the third is dynamical because it 
involves the Hamiltonian H. In terms of C there is a formal construction of an operator e 2 : 

e Q = CV, (2) 

where the operator Q is Hermitian in the conventional Dirac sense 3[. Using the operator 
Q it is possible (at least in principle) to map the Hamiltonian if to a Dirac-Hermitian 

n a 

Hamiltonian H by means of a similarity transformation of the form [8|, [9( : 

H = e~ Q l 2 He Q ' 2 . (3) 

While the similarity transformation ([3]) is isospectral, it is not unitary because Q is Her- 
mitian rather than anti-Hermitian. Hence, while H and H have the same energy eigenvalues 



E n , their energy eigenkets are not unitarily equivalent. Rather, the eigenkets \n) and \n) of 
H and H satisfy 

H\ii) = E n \n) (4) 

and 

H\n)=E n \n), (5) 

and are related by the mapping 

|n) = e e/2 |fi). (6) 

The energy eigenbra states corresponding to these eigenkets cannot be obtained from 
the kets by simple Dirac conjugation. To construct the eigenbra states of H we take the 
Dirac-Hermitian conjugate of (jlj): 

(n\H = E n (h\. (7) 

We then define 

(n\ = (n\e Q/2 (8) 

and note that (n\ is not an eigenbra state of H. Rather, the eigenbra state of H is given by 

(n\e- Q H = (n\e' Q E n . (9) 

The eigenbra and eigenket states of the Dirac-Hermitian Hamiltonian H obey the usual 
statements of orthogonality, completeness, and Hamiltonian operator reconstruction: 

(n\m) = 6 m , n , (10) 

]T|n)<n| = l, (11) 

n 

H = Y,\n)E n {h\. (12) 

n 

The similarity transformations in ([6]) and (jSJ) imply that for the non-Hermitian Hamiltonian 
H the corresponding statements are 

(n\e' Q \m) = 6 m , n , (13) 

^| n )(n|e- Q = 1, (14) 

n 

H = J2\n)E n (n\e- Q . (15) 



The norm in (TT3|) is relevant for the "PT-symmetric Hamiltonian H, with (n\e~ Q rather 
than (n\ being the appropriate energy eigenbra. In any VT theory for which the operator 
e~ s exists, there will be a positive norm of the form given in (1131) and no states of negative 
norm (ghost states). Furthermore, because [H,CVT] = 0, the PT-symmetric Hamiltonian 
H will generate unitary time evolution even though it is not Hermitian. 

To underscore the need for a non-Dirac norm for non-Hermitian Hamiltonians, we recall 
the connection between the Schrodinger and Heisenberg representations. Specifically, in 
the Schrodinger representation one introduces time-dependent states and time-independent 
operators. Ordinarily one does this for a Hermitian Hamiltonian, and it is instructive to see 
how things change in the non-Hermitian case. We thus consider Schrodinger equations for 
ket and bra states: 

•^|ob(*)> = H\a s (t)), ~ij t (Mt)\ = (Mt)\rf- (16) 

For time-independent H and H t the solutions to these equations are 

M*)) = e~ iHt \a s (0)), («s(t)| = (u s (0)\e lHH . (17) 

We introduce a time-independent (Schrodinger) operator A$ with matrix element 

(a s (t)\A s \a s (t)) = (a s (0)\e iHh A s e- iHt \a s (0)), (18) 

and define the time-dependent (Heisenberg) operator 

A H (t) = e lHH A s e- iHt . (19) 

The operator A H (t) obeys 

ij t A H (t) = A n (t)H-H^A H (t). (20) 

Since H and H ' are different when H is not Dirac Hermitian, the time derivative of An(t) 
is not given by the commutator of An(t) with H . However, regardless of whether or not 
the Hamiltonian is Dirac Hermitian, the Hamiltonian is the generator of time translations. 
Thus, in the Heisenberg representation the time derivative of an operator is always given by 

ij t A n (t) = A n (t)H-HA n (t). (21) 



To construct the Schrodinger states we start from (I2TI) and work backward. Thus, we 
replace (fT6l) and ffTTl) by 

^ |as(t)) = H \ as ^> _i dt^ s(t)l = ( d s(t)|i?, (22) 

so 

\ as (t)) = e-* m \a s (0)), (a s (t)\ = (a s (0)\e iHt , (23) 

where the bra state (as(t)\ is not the Dirac conjugate of the ket state \as(t)). In ff22l) we see 
that H acts to the right on |as(^)) and to the left on (dg(t)|. The appropriate inner product 
is given by (/3s(t)\as(t)) . Because of ( 1231 we have 

S (t)\a s (t)) = 0s(O)\e lHt e-* Ht \a s (O)) = (/%(0)|as(0)), (24) 

which is the statement of unitary time development. By contrast, the Dirac inner product 

<Ps(t)\aa{t)) = (p s (0)\e iHh e- iHt \a s (0)) ± (/%(0)|os(0)> (25) 

is not time independent. Thus, for non-Hermitian Hamiltonians unitarity is achieved by 
using an inner product that is different from the usual Dirac inner product. For VT- 
symmetric Hamiltonians the states (j3s{t)\ an d (Ps(t)\ are related by (fis{t)\ = (/3s(t)\e~ Q . 

Several models have been discussed in the literature for which one can calculate the 
operator Q in closed form. Amongst them is the Lee model [10J], where by calculating the 



correct inner product one can show that the model is explicitly ghost free ll|] . Examples in 



which one can find exact expressions for the equivalent Dirac-Hermitian Hamiltonian H in ([3 



associated with a given PT-symmetric Hamiltonian may be found in Refs. 12|, [l3j, [l4j, Il5l.ll6| 



and in our work in Ref. 5J. 

We emphasize that if we are given a Dirac-Hermitian Hamiltonian H and we convert it to 
non-Hermitian form H by means of the similarity transformation H = e Q ^ 2 He~ Q ^ 2 , we know 
that we are dealing with a Hermitian Hamiltonian in disguise. However, if we start with 
a non-Hermitian Hamiltonian H, we do not immediately know if H is a disguised Dirac- 
Hermitian Hamiltonian, and the advantage of VT symmetry is that it provides a diagnostic 
for determining whether this might in fact be the case. It may happen that the operator 
Q simply does not exist, and when this is the case, the VT- symmetric Hamiltonian will 
have no Hermitian counterpart. This situation arises when H has some complex eigenvalues 
(its VT symmetry is broken), and thus there is obviously no Hermitian counterpart. A 



transition from an unbroken to a broken VT symmetry has actually been observed in recent 
laboratory optics experiments [17| . 

However, in this paper we encounter a more serious and fundamental obstacle in trying 
to construct a Hermitian Hamiltonian H associated with a non-Hermitian PT-symmetric 
Hamiltonian H . Specifically, even if the eigenvalues of H are all real, H may be a Jordan- 
)lock matrix that has fewer eigenfunctions than eigenvalues, and thus it is not diagonalizable 
181 ] . This is the case with the equal- frequency Pais-Uhlenbeck oscillator model. For this 
case the operator Q of the unequal-frequency oscillator model becomes singular in the equal- 
frequency limit. 

III. REVIEW OF THE UNEQUAL-FREQUENCY PAIS-UHLENBECK OSCILLA- 
TOR MODEL 

The Pais-Uhlenbeck oscillator was introduced in 1950 as a simple model to explore the 
structure of a quantum system whose Lagrangian depends on acceleration as well as on 
position and velocity 4j. The Pais-Uhlenbeck action is 

J PU = 1 j dt [z 2 - (lu 2 + wf ) z 2 + ujfu^z 2 ] , (26) 

where 7, Ui, and U2 are all positive constants. Because the action depends on the accelera- 
tion, the differential equation of motion 

^ + (w? + ul)^ + u\u\z = (27) 

is fourth order. 

The Pais-Uhlenbeck oscillator model is interesting because a fourth-order wave equation 
leads to a Green's function G(E) whose denominator in energy space is quartic in the energy. 
To find the Green's function, we replace the right side of (1271) by a delta-function source 
term and take the Fourier transform. The result is 

G{E) = (E? - u,*)(E* - ulY (28) 

The advantage of such a Green's function is that it leads to Feynman integrals that are more 
convergent than the corresponding integrals constructed from propagators having quadratic 
denominators. 



One may worry that there is a price to pay for such good convergence because the form 
of this Green's function seems to imply the existence of a ghost state when ui 7^ io 2 . The 
argument goes as follows: In partial fraction form the Green's function 

G ^ = ^(^-^|) (29) 

describes the propagation of two kinds of states, one of energy uj\ and the other of energy 
uj 2 . Assuming without loss of generality that oj\ > u 2 , it appears that u 2 is associated 
with a state of negative probability because its residue contribution to the propagator is 
negative. This appears to violate the positivity condition on the weight function of the 
Lehmann representation. (Recall that when the two-point Green's function is expressed in 
Lehmann-representation form the requirement that all quantum states have positive Dirac 
norm implies that the residues of all intermediate propagating states must be strictly posi- 
tive. See, for example, Ref. [H|.) In the past, theories having fourth-order wave equations 
have been abandoned because they were thought to violate the Lehmann-representation 
positivity condition. 

However, the Pais-Uhlenbeck oscillator was recently revisited [5|] and it was shown that 
there is a realization of the unequal-frequency model in which the Hilbert space actually 
contains no ghost states. Specifically, in Ref. 5] it was shown that the above Green's- 
function argument has a subtle flaw, namely the presumption that the Hamiltonian for the 



model 



20 



21 



22 1 



2 

EX Px , 7 / 2 1 2\ 2 7222 /on\ 

Hpu = — + p z x + - (uj 1 + uj 2 ) x - -u x u 2 z (30) 



Px , . 7 / 2 . 2\ 2 7 

— +p z x+-{u l +u 2 )x --c 

is Dirac Hermitian and that its associated norm is the standard Dirac norm. [Because the 



action of (T2"6T) is constrained [20|, |2l|, |22j, to construct the Hamiltonian it was necessary to 
replace z by a new and independent variable x in the action (|26|) .] It was shown in Ref. [5| 
that one should interpret the Hamiltonian Hpu as a member of the class of non-Dirac- 
Hermitian Hamiltonians that are symmetric under combined space reflection V and time 
reversal T. Thus, as explained in Sec. [TTJ it is necessary to replace the Dirac inner product 
by the inner product that we argue is appropriate for the Pais-Uhlenbeck Hamiltonian. 

When we use the norm in (fT3"|) for the unequal-frequency Pais-Uhlenbeck model, this 
model becomes ghost free and unitary. Moreover, because the norm in (1131) is not the con- 
ventional Dirac norm, the relative negative sign that appears in the Lehmann representation 
in ff29l) cannot be interpreted as the residue associated with a negative-norm (ghost) state. 

8 



Rather, we show below that this negative sign is associated with an eigenvalue of the C 
operator defined in ([1]). In the Appendix we discuss this point further. 

To explore the structure of the unequal-frequency Pais-Uhlenbeck model in detail we 
make a further partial-fraction decomposition of the G(E) propagator in (liZUl) : 

^ ' ~ 2uj 1 (uj\-ujI) \E-ut ~ E + luJ ~ 2uj 2 (uj\ - wf) \E - uj 2 ~ E + uj 2 ) ' 

In (1311) there are two pole terms having positive coefficients and two having negative co- 
efficients. Whether or not these negative coefficients are associated with negative residues 
depends on the way one performs the contour integration in the complex energy plane. For 
the conventional Feynman contour, where the positive-frequency poles lie below and the 
negative-frequency poles lie above the leai-E axis, contour integration yields the Feynman 
propagator 



GF w = -hS dEe " iEtG{E) 



0(i) ( e - iu)lt e - iuJ2t \ 0(-t) (e iuit e iuJ2tS 



_l_ 11 ' 1 . (32) 

uj\ — u\ \ 2uji 2u 2 J oof — uj 2 \ 2uji 2u 2 J 

In f l32|) G F (t) describes the forward propagation of two positive-energy particles and the 
backward propagation of two negative-energy antiparticles and both of the uvdependent 
terms have negative coefficients. 

To avoid these negative coefficients we can instead choose an unconventional contour for 
which the poles at oj\ and — uj 2 are taken to lie below the real E axis and the poles at 
— u>x and u> 2 to lie above it. For this choice, contour integration yields the unconventional 
propagator 



G UNC (t) = -— fdEe- lEt G(E) 
2ixi J 



6{t) fe~ iuJlt , e iuJ2t \ 9{-t) ( e iuJlt e ~ iW2t \ 

+ -7T~. • (22) 



uj\ — uj 2 \ 2uj\ 2uj 2 J ujI — uj 2 \ 2lo\ 2uj 2 

Now all of the coefficients are positive, but G UN (t) describes the forward propagation of 
one positive-energy particle and one negative-energy antiparticle, and also the backward 
propagation of one negative-energy antiparticle and one positive-energy particle. 

The G F (t) propagator describes a system whose multiparticle energy spectrum is bounded 
below [the energy eigenvalues are E(ni, n 2 ) = (ni + l/2)uii+(n 2 +l/2)u> 2 } but some of its poles 
have negative residues. The poles of G UNC (t) have positive residues but the energy spectrum 



is unbounded below [the energy eigenvalues are E'(rii,n 2 ) = {n\ + l/2)u>i — {n 2 + l/2)uj2] 
and thus it exhibits forward propagation of negative-energy states. The G F (t) and G UN (t) 
propagators both seem to have problems, which explains why fourth-order theories are not 
thought to be viable. 

If we evaluate the Feynman path integral to construct the propagator, we obtain G F (t) 
directly and do not obtain G UNC (t). Specifically, for the Pais-Uhlenbeck action (126]) the 
path integral 

G(zi,Zi,Zf,Zf,T)= / T>(z, z) exp < i— I dt z 2 — (uf + ulj z 2 + lo 2 uo 2 z 2 > (34) 

laken over all paths having fixed initial and final velocities can be performed analytically 



22] . Taking its deep Euclidean time limit, where e lEt — > e Et , one finds that the low-lying 



energy eigenvalues are positive: E = u\ and E = uo 2 . Also, their excitations are just as 
required in the conventional Feynman contour prescription [22]. One should not expect a 
Feynman path integral to give an energy spectrum having negative energies because e~ Er 



would not be finite at large r, so we do not consider G UNC (t) further [23|]. Our task then is 
to find a physically acceptable quantum-mechanical interpretation for G F (t). 

To do this, we consider the Schrodinger eigenvalue problem associated with the Hamil- 
tonian H-p\j in ( |30|) . We set p z = —id z , p x = —id x and obtain 

1 d2 ■ d , 7/ 2 , 2n 2 7 2 2 2 



IX— V — (U-, + UJ 9 )X — — UJiUJnZ 



ip n (z,x) = E n ip n (z } x) } (35) 



2 7 & 2 ™dz ' 2^^ 2 ^ 2^ 2A 
which has ground-state energy E = (u\ + u; 2 )/2 and corresponding eigenfunction 



ipo{z,x) = exp 



-{oji + uj 2 )ujxU2Z 2 + i^uj^zx - -(uji + uj 2 )x 2 



(36) 



As z — > ±oo, this eigenfunction diverges. Thus, in the space of such an eigenfunction the 
operator p z = —id z cannot be Dirac Hermitian. For a canonical commutator of the form 
[z,p z ] = i, one can only associate p z with a differential operator — id z when the commutator 
acts on test functions i/j(z) that are well behaved at large z. Since this is not the case for 
the eigenfunction ipo(z, x), we see that the operator p z is not Hermitian on the real-;? axis. 
The eigenfunction ipo(z, x) vanishes exponentially rapidly for large \z\ when z is imaginary 
(or, more generally, when z is confined to the two Stokes wedges |Im(^)| > |Re(^)|). To 
exploit this fact, we perform an operator similarity transform of the quantum- mechanical 
operators z and p z : 

y = e *P**/2 ze -*P**/* = _ iz% q = e *Pzz/2 pze -*p z z/2 = ip ^ ^ 

10 



The commutator of the operators y and q still has the canonical form [y, q] = i. In terms of 
the operators y and q the Hamiltonian now takes the form 

2 

H = j- - iqx + 1 (uj\ + u 2 2 ) x 2 + |^ 2 V , (38) 

where for notational simplicity we have replaced p x by p. 

Since H and ifpu are related by a similarity transform, they both describe the same 
physics, and one can use the Hamiltonian H to explore the structure of the Pais-Uhlenbeck 
model 24J. In (138]) the operators p, x, q, and y are now formally Hermitian on the real- a; 
and real-?/ axes, but because of the — iqx term H has become complex and is manifestly 
not Dirac Hermitian. This non-Hermiticity property was not apparent in the original form 
of the Hamiltonian ifpu given in (130]) . and it is the key to uncovering the structure of the 
unequal-frequency Pais-Uhlenbeck model. 

While H is not Hermitian, with the V and T quantum-number assignments 



(39) 



we see that H is PT symmetric. Thus, H can be transformed to a Hermitian Hamiltonian 
by means of the similarity transformation H = e~ Q ' 2 He Q ' 2 . In Ref. [5] the operator Q was 
calculated exactly: 

Q = apq + j3xy, a = log 

It leads to 





p x q y 


V 


+ + 


r 


- + + - 


VT 


+ - + - 



7o>i u 2 



o>i + u 2 
UJ\ —uj 2 



j3 = aj 2 uj 2 uj- 



2, ,2, ,2 
2 



(40) 



y = y cosh 8 + i(a/(3) 1 ' 2 p smb. 8, p = pcosh 8 — i(/3 / a) 1 ' 2 ysinh 8, 
x = xcosh 8 + i(a//3) 1 ' 2 gsinh 8, q — gcosh 8 — i({3 / a) 1 ' 2 xsrnh 8, 

fr P 2 (f 7 2 2,7222 

H ^—2 + -uj 1 x+-u 1 uj 2 y, 



(41) 



2 7 2^u( 

where 8 = (a/3) 1/,2 /2, y = e _2//2 ye 2//2 , etc. The transformed Hamiltonian if is a manifestly 
positive-definite operator. Equations ( ITBl - ( TT5l) now follow directly and we obtain a viable 
unequal-frequency theory [25] . 
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IV. FOCK SPACE IN THE UNEQUAL-FREQUENCY CASE 

To construct the Fock space associated with the unequal-frequency Pais-Uhlenbeck model, 
we construct the Heisenberg equations of motion associated with H in 



y = -ix, x = -, p = iq- 7(0;? + wf )z, q = -^u)\u\y. (42) 

7 

All four operators y, x, p, and q obey the same fourth-order operator differential equations, 

"• V 1 / 2 1 2\® V 1 22 n ® X , 2 , 2\" J x , 22 n 

— + K + ^2)^0- + ^1^22/ = 0, -^1 + K + ^2)^2- + ^1^2^ = 0, 

"P \ ( 2 1 2\^ P 1 22 n " 5 1 / 2 1 2\ « <7 , 2 2 n //io\ 

— + K+^ 2 )— +w 1 w 2 p=0, — + K+^ 2 )— + Wl w 2 g = 0, (43) 

which are identical in form to the classical Pais-Uhlenbeck equation of motion. 

Despite the presence of the complex number % in H, all the coefficients in the quantum- 
mechanical equations of motion are real. From (1431) each of the dynamical operators can be 
expanded in terms of a basis consisting of two raising and two lowering operators having 
e ±iuj!t anc ^ e ±iu 2 t tj me dependence. From (142 p we thus obtain 



y = -ia ie - iuJlt + a 2 e~ iuJ2t - ia ie iuJlt + a 2 e iu>2 \ 

x = -iujxaxe-^ 11 + u 2 a 2 e- iuJ2t + itu^e^ - cu 2 a 2 e iuJ2 \ 

p = 7 [_o; 2 aie -^i* _ iu 2 2 a 2 e- iuJ2t - ujla ie iullt - iu 2 2 a 2 e iuJ2 % 

q = 7ajicj 2 [-cj 2 aie^ lt - itu^e^ 2 * + tu 2 a ie iuJlt + zuw^ 2 *]. (44) 

The operators di and d 2 are not the Dirac-Hermitian adjoints of the lowering operators Oi 
and a 2 because they are the raising operators for a non-Hermitian Hamiltonian. In terms of 
the raising and lowering operators, H takes the diagonal form 

H = 2 r y{u\ - u 2 ) (ufaiai + uj\a 2 a 2 } + -{u\ + u 2 ), (45) 

where the operator commutation algebra is given by 

[a x , oi] = [2ju>!(ujI - ul)]' 1 , [o 2 , d 2 ] = [2jlj 2 (u}1 - u^)] -1 , 

[a 1; a 2 ] = 0, [oi, d 2 ] = 0, [di, o 2 ] = 0, [di, d 2 ] = 0. (46) 

In ( 1451) and ( 1461) the relative signs are all positive, so these equations define a standard 



two-dimensional harmonic oscillator system 26J]. The Heisenberg equations of motion for 
the raising and lowering operators are 

di = i[H,ai] = —iuidi, di = iuiai, a 2 = —iuo 2 a 2 , d 2 = iuo 2 a 2 . (47) 

12 



The solutions to these equations are Oi(t) = Oi(0)e" 



■iui\t 



6i(t) = ai(0)e iwit , a 2 {t) 



a2(0)e tuJ2t , a 2 (t) = a 2 (0)e tuJ2t , just as is required of any set of raising and lowering operators. 

V. FOCK SPACE IN THE EQUAL-FREQUENCY CASE 

The apparatus of VT quantum mechanics, as developed above, is readily implementable 
for the unequal-frequency Pais-Uhlenbeck model, but the operator Q, as given in PHI) , 
becomes singular in the equal- frequency limit uj\ — uj 2 —>■ 0. Additionally, in this limit the 
partial-fraction form of G(E) in (1291) is not valid. Thus, the equal-frequency limit of the 
Pais-Uhlenbeck model must be treated separately from the unequal-frequency case. We 
will see that because Q becomes singular, the PT-symmetric Hamiltonian H develops a 
nondiagonalizable Jordan-block structure 27|. Because this happens it is no longer possible 
to construct an equivalent Hermitian Hamiltonian H. Moreover, the VT sector of the 
equal-frequency theory becomes an independent and self-contained realization of quantum 
mechanics. An objective of this paper is to show that in the singular equal-frequency limit 
the unitarity of the Pais-Uhlenbeck theory is not lost. 

To obtain the equal-frequency limit of the Pais-Uhlenbeck theory we must find a basis in 
which, unlike (1451) and (1461) . the operators and operator algebra are continuous in the limit. 
To do this we define 

uj\ = uj + e, u 2 = 00 — e 



(48) 



and introduce the new operators 

e 



a 



01 1 + 



01 1 + 



2a; 
e 



+ ia 2 11 - 
+ ia 2 [l — 



e 

e 

2u 



2a; 

€ 

2a7 



ai — ia 2 ) 



ai — ia 2 . 



(49) 



These new operators obey the commutation algebra 



[a, a] = A, [a, b] = //, [b, a] = /i, [b, b] = A, [a, b] = 0, [a, b] = 0, 



(50) 



where 



A 



2a? 



/' 



167(a; 2 — e 2 )a; 3 ' ' 167(a; 2 — e 2 )a; 3 
In terms of these new operators the unequal- frequency position operator y(t) is 



y^o 



-iujt 



-i(a — b)cos et 



2bu 



-sin et 



+ e 



■iuit 



-i[a — b)cos et -\ sin et 



(51) 



(52) 
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while the unequal-frequency Hamiltonian H in (|38|) is rewritten as 

H € ^ = 8<yu 2 e 2 (aa - lb) + 87CJ 4 (2bb + ab + So) + w. 



(53) 



The advantage of these new operators is that -ff(e) and y(e), and also p(e), x(e), and g(e), 
are continuous in the e — ^ limit: 



(54) 

y e =o = e "~ [— i{a — 0) — 'iouii\ -+- e"~ — 2^a — oj + zoa>t . (55) 

The commutation relations (j50p are also continuous in the limit e — > and, together with 
the commutation relations with H(e), they tend to the gaugelike form 



H e=0 = 87^(266 + ab + ba) + u, 
' iut \-i(a -b)- 2but] + e iu}t \-i(a -b) + 2but 



\a,a\ 



0, [a, 6] 



:, [b,a] 



57CJ Q'JW 

[H £=0 , a] = o>(d + 26), [H e=0 , a] = -uj(a + 26), 



[6, 6] = 0, [a, 6] = 0, [a, 6] = 0, 



[H e=0 ,b] = tub, [H e=0 ,b] 



-uib. 



(56) 



Note that when 6 ^ the states a\Q) and b\Q), where \Q) is the no-particle vacuum 
annihilated by a and 6, are not eigenstates of H e ^Q. Rather, the action of the unequal- 
frequency Hamiltonian on these states is given by 



H e ^a\Q) = — \(4lu 2 + e 2 )a\n) + (Alo 2 - e 2 )6|ft) 

H e ^ l\Q) = — [eV|fi) + {4u 2 - e 2 )b\fl) 
2ui l 



and the Hamiltonian acts in the one-particle sector as the nondiagonal matrix 

,2 , r 1 Al ,2_Jl\ 



M, 



e^O 



I /4a; 2 + e Auj 2 - r 



2U V e 2 4^2 _ £ 2 J 

In the one-particle sector the Hamiltonian has two eigenstates 

H e ^ \2uj ± e) = (2a; ± e)\2u ± e), 

where 

\2u ± e) = [±eo + (2cj t e)b] |fi). 

A similarity transformation diagonalizes M e ^ : 

'4a; 2 + e 2 4cj 2 - e 2 



s- 1 



2u 



2 ^2 



5 



Au z -e 



2u + e ' 
2a; -e 



(57) 



(58) 



(59) 



(60) 



(61) 
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where 

1 (2uj + e -(4cj 2 -e 2 )e~ 



S 



26^/2(20; + e)V2 ^ e ( 2w + e)e 2 . 

1 ({2u + e)e 2 (4cu 2 -e 2 )e' 



5_1 2eo;^(2a; + e)V* 1 - e 2u; + e J ' ^ 

Now let us examine the limit e — > 0. From (16"Uj) we can see that as e — > the two 
eigenstates \2u ± e) collapse onto one state b\Q). Thus, -H"(e) loses a one-particle eigenstate 
in this limit. This is to be expected because the commutator [H e= o, a] contains components 
parallel to both a and b. To understand why an eigenstate has disappeared, we note that as 
e — > the matrix M(e) in (I58p takes the upper-triangular Jordan-block form 

M £=0 = 2w f V (63) 

The matrix M £=0 possesses two eigenvalues, both equal to 2u, but it has only one eigenstate 
because the eigenvector condition 

c :)o-m-o 

permits only one solution, namely that with d — 0. Despite not being Dirac Hermitian and 
despite having lost an eigenstate, H e=Q continues to be VT symmetric under the transfor- 
mations (1551) . and thus its eigenvalues remain real. Furthermore, as e — > 0, S and S 1-1 in 
(|62|) both become singular, so it is not possible to diagonalize M e= Q. 

This same pattern occurs for the higher excited states of H. (For example, one can readily 
show that the three two-particle states collapse onto one common eigenstate as e — > 0.) To 
discuss the e — > limit for arbitrary multiparticle states, we return to the operator Q in 
(|40p and note that at t = we can use (|44p to write 

pq + ~i 2 iJ\ij)\xy 

= ry 2 co>iCo> 2 (^ 2 - ojD [{u)\ - uj 2 )(a 1 a 2 - Oi0 2 ) + (^i + uj 2 )(a 1 a 2 - o 2 Oi)] . (65) 

In the e — * limit we thus obtain 

pq + ^ujfujjxy -> 87 2 w 5 (6 2 -b 2 + ba- 66). (66) 

The quantity pg + ^ufu^xy is well-behaved as e — > 0, and the singularity acquired by the 
coefficient a in (1401) cannot be canceled. The similarity transform e _S//2 in ([3]) also becomes 
singular. 
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Since the operator Q creates two-particle pairs in (jfi"5|) and (I6"6"j) . e _S//2 is a singular 
operator in every multiparticle sector of H, with each such sector developing Jordan-block 
structure in the limit. The equal-frequency Pais-Uhlenbeck model possesses no equivalent 
Dirac-Hermitian counterpart; as we show in the next sections, it constitutes a self-contained 
realization of quantum mechanics that exists in its own right. 



VI. EIGENFUNCTIONS OF THE UNEQUAL-FREQUENCY THEORY 



The outstanding property of a Hamiltonian in Jordan-block form is that its eigenstates 
are incomplete. One may ask, where do the other eigenstates go in the equal-frequency 
limit, and how can one formulate quantum mechanics if the space of energy eigenstates is 
incomplete? To answer these questions for the Pais-Uhlenbeck model, we must construct the 
eigenfunctions of the unequal-frequency theory, where there are no Jordan-block structures, 
and we must then track what happens to the eigenfunctions in the equal- frequency limit. 

Noting that the coordinate-space representation of the Hamiltonian ( l38l) is not symmet- 
ric because it has a term that acts like xd y , we will need to distinguish between a right 
Hamiltonian, which operates to the right, and a left Hamiltonian, which operates to the 
left. When acting to the right on well-behaved states, the commutator [y, q] = i is realized 
by setting q = —id y , but when acting to the left it is necessary to use q = +id y . We thus 
obtain the right and left Schrodinger equations 



.dtjj^(x,y,t) 
1 dt 

■ d^n(x,y,t) 
dt 



27 dx 2 
1 d 2 



O 1 7/ 2 1 2\ 2,7222 

~d~ 2 2 UlU2y 



^n( X ,V^), 



(67) 



d 



/ / 2 2\ 2 / 2 2 2 



The ground state, whose energy is E 



IW! 



&(x,y,t). (68) 

co> 2 )/2, has right and left eigenfunctions 



*l>o(?,y,t) 



exp 



7, 



7 , 



—h{u\ + u 2 )uJiuj2y 2 - -fuJi^yx - ~(u>i + u 2 )x 2 - ~(Ul + u 2 )t 



4>o(x,y,t) = exp 



7, 



7, 



-~{u\ + u) 2 )oji uj 2 y 2 + ^uj!UJ 2 yx - -(wi + u 2 )x 2 + ~(Ul + UJ 2 )t 



(69) 
(70) 



For this ground state we can define a normalization integral of the form 



iV n 



dxdyipv(x } y } i)il%(x,y } t) 

dxdyexp — 7(0*1 + u^u^^ 2 — 7(0*1 + u) 2 )x 2 



TV 



7(01 + Co>2)(o , lO> 2 ) 1 / 2 



(71) 
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with 7Vo being time independent, finite, and real. 

The two one-particle states with energies E\ = E + uji, E 2 = E + u 2 have eigenfunctions 

i>?{x, y, t) = (x + u 2 y)^(x, y, t^ 1 ^, ^(x, y, t) = (x - cj 2 y)^(x, y, t)e^, 
ijf(x, y, t) = (x + u iy )^(x, y, t)e~^\ $(x, y, t) = (x - u iy )$(x, y, t)e^\ (72) 



whose normalization integrals are 

7r(cJi — UJ 2 ) 



N- L 



3/2 1/2 ' 



N, 



n(ui — uj 2 ) 



(73) 



2^ 2 (uoi + uj 2 ) 2 lUi /z uj2 Z Z^i^i + 0O 2 ) 2 UJi UJ 2 

Again, these normalizations are time independent, finite, and real, but N 2 is negative, a 
crucial issue that we return to and resolve below. 

Analogously, the three two-particle states with energies E 3 = E + 2ui, £4 = E q + uji + uj 2 , 
and E$ = Eq + 2uj 2 have right and left eigenfunctions 

1 



and 



*l>f(x,y,t) -- 
i^{x,y,t) 



(x + u 2 y) - 



2^fooi 
(x + uxy)(x + u 2 y) 

(x + wij/) 2 



^(x,y,t)e- 2i ^ 
1 



2juj 2 



7(^1 + L0 2 ) 

^ R (x, 2 /,t)e- 2 ^ t 



^(x.y.Oe-***-***, 



(74) 



(x - u 2 yf - 



1 



270;! 
(x - uj x y)(x - u 2 y) - 
1 



^{x,y,t)e 2i ^ 
1 



(x - u^y 



2^oj 2 



7(cji + u 2 ) 
^(x,y,t)e 2 ^. 



i}> (x,y,t)e 



iu>it+iu>2t 



(75) 



The normalization integrals are time independent: 

7r(cJi - U 2 ) 2 



N 3 

N 5 



2^^ + ^{u^y/ 2 ^ 
w(uji — uj 2 Y 
47 3 (cji + uo 2 Y(uo 1 uo 2 ) 1 l 2 uo 1 uo 2 ' 
7t(cji - UJ 2 f 



(76) 



27 3 (cji + cj 2 ) 3 (ccW 2 ) 1/2 ^2' 
Note that JV4 is negative. 

To relate these normalization integrals to the Hilbert-space norms of the eigenstates of 
H, we recall that the energy eigenstates of the Hamiltonian are eigenstates of both the 
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VT operator and of the C = e Q V operator introduced in Sec. [TT1 The general procedure 
for constructing such states for symmetric Hamiltonians is given in Ref. |3y, and we adapt 
it here for the nonsymmetric case. Since the VT operator is antilinear, in general its 
eigenstates have eigenvalue e ia , where a is a real phase that depends on the eigenstate. 
Multiplying each eigenfunction by e ia ^ 2 gives a new eigenfunction that is still an eigenstate 
of the Hamiltonian, but now it has VT eigenvalue equal to one. Because the operator C 
obeys the algebraic conditions in ([I]), its eigenstates are also eigenstates of VT and of H 
and all of its eigenvalues C n are ±1. Consequently, each energy eigenstate is an eigenstate 
of CVT with eigenvalue C n = ±1. 

When the Hamiltonian H is symmetric, the completeness and normalization conditions 
have the form 

J2 e^CVT^x', y', t = 0)]V«(s, V, t = 0)e~^ 4 = 8(x - x')5{y - y'), 

n 

Jdxdye iE " t [CVT^ n (x,y,t = 0)}ip m {x,y,t = 0) e -^ f = 8 m>n , (77) 

where the summation is taken over all the eigenstates of the Hamiltonian. For the Pais- 
Uhlenbeck oscillator, H is not symmetric and thus we must distinguish between right and 
left wave functions. We identify ipn(x', y', t = 0) = CVTip^(x', y', t = 0) and then assign the 
C n eigenvalues as follows: We take Cq = 1 for the ground state ipo(x,y,t), C\ = 1 for the 
one-particle state ipi{x,y,t) and C 2 = —1 for the one-particle state fa(x, y, t). The three 
two-particle states in (JTlj) then acquire eigenvalues C3 = 1, C4 = — 1, C5 — 1. 

The alternations in sign of the C n parallel the alternations in sign of the normalization 
integrals given in ( I7T1) . (|73l) . and ( !76|) . with the same pattern repeating for the higher 
multiparticle states; that is, negative signs occur when the number of 0J2 quanta is odd. 
Thus, for any eigenstate the sign of C n is precisely the same as that of its normalization 
integral N n = J dxdy ipn(x, y, t)ip^(x, y, t). Consequently, from (1771) the correct completeness 
relation and normalization conditions for the states of the theory are 

J2 tt(x', y', t)£-r1>*{x, y, t) = 8(x - x')S(y - y'), (78) 

/ dxdyipn(x,y,t)—^-ril} m (x,y,t) = 5 m , n . (79) 

Both this CVT norm and the Fock space norm (n\e~ Q \m) = 5 m ^ n in (j!3p are positive. Thus, 
there are no negative-norm states in the unequal-frequency Fock space and the unequal- 
frequency Pais-Uhlenbeck model is unitary. 
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Because the unequal-frequency energy eigenstates are complete, we can expand an ar- 
bitrary wave function in terms of them as ip R (x,y,t) = J2 n a n ip R (x,y,t) and ip L (x,y,t) = 
J2n a n' l l ! n{ x iyit)- Then, since the energy eigenstates form an orthonormal basis with real 
energy eigenvalues, it follows that the norm / dx dy tp L (x, y, t)Ctp R (x, y, t) is both time inde- 
pendent and real. Therefore, probability is preserved and the unequal-frequency theory is 
unitary under time evolution. 

Given the norm in (1791) and writing the norm in ffT3|) in the form (n\e~ Q \m) = 
J dxdy (n\e~ Q \x,y){x,y\m), we make the identifications 

C 1 

I AH 1/2 ^"fo y ' *) = ( n l e l x » V, t), . u/2 ^(^ v, t) = (s, 2/> t\n). (80) 

Then, by using the relations in (IHUI) and substituting i7 = J2 n \ n )E n {n\e~ Q into the expres- 
sion G(x, y, x', y', t) = (x, y,t = 0\e~ tHt \x', y', t = 0) for the propagator, we obtain 

G(x, y, x', y', t) = £ ift(x, y, t = 0)J^- e - iE ^(x', y', t = 0). (81) 

Thus, the negative sign of C2 accounts for the negative sign of the c^-dependent term in 
the propagators of (129]) and (1321) . This shows that the good convergence associated with 
fourth-order propagators can be achieved in a Hilbert space whose norms are all positive 
and need not be in conflict with the requirement of unitarity. 

VII. EIGENFUNCTIONS OF THE EQUAL-FREQUENCY THEORY 

In the equal-frequency limit the Hamiltonian (138]) takes the form 



2 
H = — -iqx + 7U>V + -wV • (82) 



Its ground state has energy E = u, left and right eigenf unctions 

Vfr(x, y, t) = exp — *yuj 3 y 2 — ^uP'yx — ^ujx 2 — iut 



ifj Q (x, y, t) = exp —"fuoy +70; yx — ^ux + iut , (83) 



and a normalization integral 



7[ 



iV = / dxdyt(jo(x,y,t)ij R (x,y,t) = - — -, (84) 



which is time independent, finite, and real. 
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The equal-frequency theory differs from the unequal-frequency theory in that there is 
only a single one-particle eigenstate instead of two one-particle eigenstates. The energy of 
this state is E\ = 2u and its eigenfunction is 

4>?(x, y, t) = {x + uy)$(x, y, t)e~^\ $(x, y, t) = (x - uy)^{x, y, t)e iu)t . (85) 

In the equal-frequency limit the two unequal-frequency eigenstates ipf(x, y, t) and i$(x, y, t) 
of (1T2"j) collapse onto one state ipf(x, y, t) = [ipf(x, y, t) + ipf(x, y, t)]/2 (and likewise for the 
left eigenfunction). 

The disappearance of eigenstates in the equal-frequency limit is generic. This same 
collapse of eigenstates occurs for the higher excited eigenstates, with the three unequal- 
frequency two-particle eigenf unctions in ( f74l) and (1751) collapsing onto a single equal- 
frequency second excited eigenstate having energy E 2 = 3u: 

1 



4>f(x,y,t) 



(x + ujy) - 



27a; 



'.R/„ „. +\„-2iuit 



(x - u)yf 



ip (x,y,t)e 

^(x,y,t)e 2i ^. (86) 



27a; 

Even though the normalization integral iVo for the equal-frequency ground state is posi- 
tive, for the equal-frequency first excited state we find that 

Ni = fdxdyip];{x,y,t)i>?{x,y,t) = fdxdy(x 2 - w y) e -27"V-2 7 ^ 2 = . (87) 



The norm of this state vanishes because the unequal-frequency normalization integrals N\ 
and N 2 in (1731) both vanish in the equal-frequency limit, and the unequal-frequency states of 
(I72|) are orthogonal before the limit is taken. Thus, tpf(x,y,t) ends up being orthogonal to 
itself, that is, to ^\(x,y,t). This same situation repeats for the higher excited states with 
all of the two-particle norms of ( ITS"]) collapsing onto ^2 = 0. 

The emergence of zero-norm states is characteristic of Hamiltonians having Jordan-block 
structure. For the two-dimensional Jordan-block matrix (I64p . the right eigenvector is given 
by the column (1,0), while the left eigenvector is given by the row (0, 1) rather than by 
the row (1,0). Thus, the norm, which is the product of the left and right eigenvectors, is 
zero. Yet, despite the presence of zero norms, we will show below that probabilities in the 
equal-frequency theory are still nonzero. 

The energy eigenstates that disappear in the equal-frequency limit are replaced by an 
equal number of nonstationary states. To demonstrate this phenomenon, we form linear 
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combinations of the unequal-frequency eigenf unctions ipf ( x , V, t) and ipf(x, y, t) with coef- 
ficients that behave like 1/(^1 — 0J2) and then track the limit. To do this, we extract the 
terms in ipf(x, y, t) and 4>f(x, y, t) that are linear in uj\ — u 2 . This yields the nonstationary 
equal-frequency state 

fc{x, y, t) = lim ^iM__jgfeM = [{x + uy)lt + y] ffifa y , t)e -^ 
fa{x,y,t) = Km ^ {x ' V ' t) -rt {x ' y ' t) = [-{x-uy)it-y]^{x,y,t)e^. (88) 

Since this state is not stationary, it is not an eigenstate of H. However, it does satisfy the 
time- dependent equal-frequency Schrodinger equation: 

z-^ R (x, y ,t)=(-^—-x— + 7 wV + |u,yj ^ R (x, y, t), 

-^ L (*, v, t) = (--^ + *§y- + W + |«v) ^fo y, *)■ (89) 

Similarly, on expanding the unequal-frequency two-particle states of (!74|) to order {u\ — cu 2 ) 2 , 
one obtains [in addition to i>^(x, y, t)] the two states 



e^o 2e 



(x + uyf - - — ) 2it + 2xy + 2cuy 2 - 



27a;/ 270;^ 



^(x,y,t)e-^, (90) 



and 



?r/ ,x r 2^ R (x,y,t) - ipf (x,y,t) -i/jf(x,y } t) 
^2b{x,y,t) = km — 

(x + uyf -- — ] 2t 2 - \2xy + 2ujy 2 - - — -) lit - 2y 2 + 



27a; J \ 2ju 2 J 27a; 3 

x^ R (x,i/,t)e- 2 ^. (91) 

Both of these states satisfy the equal-frequency time-dependent Schrodinger equation even 
though neither state is stationary 

The picture is now clear. The unequal-frequency states ipf(x, y, t) and ipf(x, y, t) are en- 
ergy eigenstates, but linear combinations of them are not because the states are not degen- 
erate. However, such linear combinations are still solutions to the unequal-frequency time- 
dependent Schrodinger equation. When e — ► 0, we obtain two new wave functions ipf"(x, y, t) 
and ?pf a (x, y,t), which are solutions to the equal-frequency time-dependent Schrodinger 
equation. One state is stationary, so it also solves the equal-frequency time-independent 
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Schrodinger equation and thus it is an energy eigenstate, while the other state is not. The 
same pattern repeats for the higher excited states. The equal-frequency Hamiltonian is 
a Jordan-block matrix that has fewer eigenstates than eigenvalues, and the form of the 
unequal-frequency Hamiltonian becomes Jordan block as e — > 0. The counting of wave func- 
tions is continuous in the limit, and it is just the counting of stationary states that is not 



Or* ~ r 3 r a ~ i 

i-KTj dxdyi/jB(x,y,t)^(x,y,t) = -J dxdyx — \j^(x,y,t)^(x,y,t)\ . (92) 



continuous. The nonstationary states replace the missing stationary states in the limit 29]. 

Instead of counting energy eigenstates, if we count Fock states, we see that the unequal- 
frequency Fock space built from the raising operators a\ and 02 in (|46|) and the equal- 
frequency Fock space built from the raising operators a and b in ( |56l) both have the dimen- 
sionality of a two-dimensional harmonic oscillator. The eigenspectrum of the equal- frequency 
Hamiltonian H e=0 is like that of a one-dimensional harmonic oscillator. However, the count- 
ing of states in the full Fock space is continuous in the limit e — > 0. 

Given the form of the Schrodinger equation in (1591 . any pair of its solutions ip^{x, y, t) 
and ?Pb(x, y, t) obey 

d r , , ?t ., .^ n , .x r , , d_ 

dy 

Because the product ipQ(x,y,t)ilj^-(x,y,t) = exp(— 27a; 3 ?/ 2 — 2 r yux 2 ) is exponentially sup- 
pressed, the overlap integrals for any pair of multiparticle eigenfunctions are time indepen- 
dent: 

i— J ' dxdyife(x,y,t)$i(x,y,t) = 0, (93) 

where the indices A and B are 0, 1, la, 2, 2a, 26, and so on. (These eigenfunctions have the 
form of polynomials in x and y multiplied by the ground-state eigenfunction.) Because H 
is not Dirac Hermitian, (1931 has the generic time-independent form in (T24|) . 
For modes that are energy eigenstates, (193)) reduces to 

{E A -E B )fdxdyfe{x,y,t)4%(x,y,t) = 0, (94) 

so modes having unequal-energy eigenvalues are orthogonal as usual. When A is an eigen- 
mode, the diagonal integrals / dxdyip^(x,y,t)ipj^(x,y,t) all vanish except for the ground- 
state eigenmode because all eigenstates other than the ground state have zero norm. The 
overlap integrals of the nonstationary wave functions, either with each other or with the 
stationary solutions, need not vanish because (!9B1 requires that such overlaps be time in- 
dependent, and not that they vanish. However, the equal-frequency wave functions are 
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constructed as limits of linear combinations of unequal-frequency modes, so the orthogonal- 
ity of multiparticle eigenstates with different numbers of particles in the unequal-frequency 
theory translates into the orthogonality of the associated states in the equal-frequency case. 
Thus, even though not all of the equal-frequency theory states are stationary, both of the 1 
and la states are orthogonal to all of the 2, 2a, 26 states, and so on. 

Within any given multiparticle sector the overlaps need not all be zero. For the typical 
one-particle sector, for instance, we evaluate the overlap integrals directly and find that 

dx dy i)\ (x, y, t)ij>f(x, y, t) = 0, 



a. a. f a. a 

dxdyil)\(x,y,t)i^ a {x,y,t) = / dxdy^\ a {x,y,t)^ K {x,y,t) 



TX 



8 7 2 w 4 ' 

dxdy^ a (x,y,t)ipf a (x,y,t) = -r-5-5. (95) 

87^0;° 

This confirms that these overlaps are indeed time independent. The overlap integral 
f dx dyipi(x,y,t)ipf a (x,y,t), for example, is time independent because the coefficient of the 
term that is linear in t is just the zero norm J dxdytp\(x,y,t)tp R (x,y,t). 

An alternative way to derive the relations in (1951) is to take the e — > limit of the 
normalization integrals Aq and N 2 in (1751) . To order e this yields 

dxdy[^(x,y)ipf(x,y)-e^ a (x,y)^f(x,y) - eip];(x,y)ipf a (x,y)] = -7-5-7, 

-S A A A ^, -S TTf 

dx dy [^(x, y)ipf(x, y) + eip^ a (x, y)ipf(x, y) + e^(x, y)ipf a (x, y)\ = -^-3^4 > ( 96 ) 

from which the relevant relations in (|95|) follow. This procedure implies that in taking the 
e — > limit of the time-independent normalization integrals of the unequal-frequency theory, 
we get time-independent expressions, some of which are nonzero, because as (175)) and (176|) 
show, the unequal-frequency-theory normalization integrals only vanish as powers of e. 

The nonstationary wave functions are not energy eigenstates, but at t — the various 
multiparticle wave functions ipn ( x 5 V-> t) an d ^ a (. x i Vi ( n = 0? 1> 2, ..., a = a, b, ...) contain 
polynomials in x and y of degree n, with just enough freedom to construct any arbitrary 
polynomial function of x and y. Hence any initial wave function ip R (x,y,t = 0) of x and 
y can be expanded in terms of a complete basis of polynomial wave functions: ip R (x, y,t — 

0) = En ani>n( X > V) * = °) + E„,a a n,a^n,a( X > J/> * = °)" 

The left wave functions have analogous expansions, with relations such as those in (|95|) 
implying that the quantity N(ip,t = 0) = / dx dy ip L (x , y , t = 0)ilj R (x,y,t = 0) is nonzero. 
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Thus, we can use the nonstationary solutions to construct any initial state whose initial 
probability N(tp,t = 0) is nonzero despite the presence of zero-norm energy eigenstates. 

Each of the i) R {x, y, t = 0) and ipna( x i Vi^ = 0) basis wave functions is a solution to the 
time-dependent Schrodinger equation, so each has a uniquely specified time evolution. Thus, 
given (I9"5j) . we see that the probability integral / dxdyif) h (x,y,t)if; R (x,y,t) is preserved in 
time, with all the terms that involve powers of t dropping out of N(ip,t). Therefore, the 
equal-frequency theory is unitary. 

To summarize the nature of completeness in the Jordan-block case, where there are 
fewer energy eigenstates than the dimensionality of the Hamiltonian, the missing eigenstates 
are replaced by an equal number of nonstationary solutions to the Schrodinger equation. 
Together the stationary and nonstationary states form a complete basis that may be used 
to construct an initial wave packet. Because these states are complete, the normalization of 
the wave packet is preserved in time. 

Terms involving powers of t contribute to the Green's functions of the equal-frequency 
theory even though they play no role in probability integrals such as N(tp,t). We construct 
the equal-frequency Green's function G(x, y, x', y', t) as the limit of the unequal-frequency 
Green's function G(x,y,x',y',t) given in (IHTl) . In the unequal-frequency theory the repre- 
sentative one-particle contribution to (!HT|) is 

G(x,y,x',y',t)( 1 - ) 

= iP R (x,y,t = 0)^-e- iE ^(x',y',t = 0) + ^ R (x,y,t = 0)-^- e - iE ^(x', y', t = 0) 



ArJi, ,2(, ,2 _ ,2\l/2 r 

= ^ R (x, y, t = 0)rt(x>, y', t = oje-W** 47 "^ €) e -**( w + e) 

x (x' — ujy + ey')(x + uy — ey) — e t€t (u — e)(x' — u>y' — ey')(x + uy + ey) . (97) 



Despite the presence of terms that behave as 1/e, the e — > limit of fl97|) exists: 

G(x,y,x',y',t) {1) -. ^ R (x,y, t = 0)^(x',y', t = 0) e - 2 ^^-^- 

x [(1 — iu>t)(x' — u>y')(x + uy) + ujy'{x + ujy) — u(x' — ujy')y] ,(98) 



and using (1851) and (1881) . we rewrite (1981) as 



2, ,3 



yy u 



G(x, y, x', y', t) {1) = [w(x, y, t)$(o/, y', t = 0) - unf^{x, y, t'm{x', y', t = 0) 

-wj>?(x,y,t)$ a (x , i y',t = QJ\. (99) 
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Equation (1991) reveals the role played by the nonstationary states in the propagator. 
Because G(x, y, x', y', t) describes the propagation of a wave packet that is localized at (x, y) 
at an initial time t, both stationary and nonstationary wave functions are contained in the 
wave packet, and both are needed to form a complete basis with which to construct localized 
wave packets [301 ]. 

When the initial and final states are at the same time, the Green's function gives the 
normalization of the eigenstates of the position operator according to (x,y,t\x',y',t) = 
S(x — x')S(y — y'). Consequently, the e — > limit of the unequal-frequency completeness 
relation given in ( 1751) must recover this property of ( 1991) . Explicit evaluation of ( 1781) yields 

fc R (x, y, t)$(x', y', t) - u4? a (x, V, t)ti(x', y', t) - ^f(x, y, t)fyx', y', t) 



IT 

„3, ,4 



$2(x, y, t)fy{x', y', t) - uji{)f a {x, y, t)fy{x', y', t) - u$f{x, y, t)fyrf, y', t) 



IT 
3, ,6 



+ — — 



^f a (x, y, t)^ a {x', y', t) - ijf b (x, y, t)^(x', y', t) - tfjf(x, y, t)^ b (x', y', t) 

TT L J 

5(x - x')5(y - y') (100) 



+ 



at all times, just as required. Using (195|) and the orthogonality of the different multipar- 
ticle wave functions, one can verify (11001) by projecting with / dxdyip\{x ) y 1 t) 1 and so on. 
Equation fllOOj) thus generalizes the standard completeness relation to the nonstationary 
case. 

We have shown in this section how to develop a consistent quantum-mechanical theory 
given that the energy eigenstates do not form a complete basis. We did not make use of 
the fact that energy is quantized. Rather we introduced a specific canonical form for the 
commutators of the position and momentum operators, and we did this without reference to 
the structure of the Hamiltonian. Demanding that position and momentum operators such 
as x, p, y, and q be quantum operators required that we specify a Hilbert space on which 
they operate. In the basis in which the position operators are diagonal, we introduced a 
complete set of basis vectors \x,y), and in so doing we specified the Hilbert space once and 
for all. 

An alternative but equivalent prescription is to represent the position and momentum 
operators as infinite-dimensional matrices acting on a Fock space, where all the states are 
created from a no-particle state. The raising operators then generate the complete set of 
basis vectors. This procedure works whether or not the Hamiltonian commutes with the 
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Fock-space number operator. This construction does not involve the Hamiltonian, and it 
does not depend on how many energy eigenstates the Hamiltonian possesses or on which 
states in the Fock space are its eigenvectors. 

Normally, in quantum mechanics the Hamiltonian is a Hermitian operator whose energy 
eigenstates form a complete basis. This basis is in one-to-one correspondence with both 
the coordinate-space basis and the Fock-space basis. Hermitian Hamiltonians can be diag- 
onalized, so for a Hermitian Hamiltonian it is advantageous to use its eigenvectors rather 
than any other set of vectors as the basis vectors. In this paper we have shown that one 
need not have a complete basis of energy eigenstates to characterize a quantum-mechanical 
Hilbert space when the Hamiltonian is in Jordan-block form. If the Hamiltonian cannot be 
diagonalized, the coordinate-space and Fock-space bases are central and the nonstationary 
solutions to the Schrodinger equation play a role that they do not play in the case of a Her- 
mitian Hamiltonian. The lack of diagonalizability of Jordan-block Hamiltonians is not an 
impediment to the construction of a fully consistent and unitary Jordan-block quantum the- 
ory, and for such cases PT quantum mechanics represents a distinct realization of quantum 
mechanics that exists in and of itself. 

VIII. CONCLUSIONS AND COMMENTS 

The fourth-order-derivative Pais-Uhlenbeck model is rich and instructive and we have 
used it to examine many issues of contemporary concern. In our previous paper J5( we 
showed that the unequal-frequency version of the model has a consistent quantum real- 
ization in which the spectrum is real and bounded below, the Hilbert space of states is 
ghost free, and time evolution is unitary. In this paper we have examined the case of the 
equal-frequency Pais-Uhlenbeck model. We have constructed this model by performing the 
equal-frequency limit of the unequal-frequency model. This limit is singular because the 
Hamiltonian develops a Jordan-block structure and many of the eigenstates of the Hamilto- 
nian disappear. Nevertheless, we have shown that the limiting theory remains a consistent 
and unitary quantum theory. 

Our solution to the equal-frequency Pais-Uhlenbeck model stems from our work on 
the unequal-frequency model. To find a physically acceptable realization of the unequal- 
frequency model we established in Ref. |5j that the Pais-Uhlenbeck model Hamiltonian H is 



2(3 



actually not Dirac Hermitian, but is instead VT symmetric. This symmetry allowed us to 
construct a similarity transformation H = e~ 2//2 ife s / 2 , which produces the Dirac Hermitian 
Hamiltonian H having the same eigenvalues as H. We calculated the operator Q exactly 
and in closed form. The Q operator reveals the singular nature of the equal-frequency 
limit because it ceases to exist in this limit. Thus, there is no equivalent Dirac-Hermitian 
Hamiltonian for the equal-frequency theory. Nevertheless, we have shown in this paper that 
all of the eigenfunctions that disappear in the equal-frequency limit are replaced by time- 
dependent solutions to the Schrodinger equation. As a result, completeness is maintained 
and the model continues to exhibit unitary time evolution. Thus, we have shown that the 
equal-frequency Pais-Uhlenbeck model is a unitary VT quantum theory that has no equiv- 
alent Hermitian counterpart, and so VT quantum mechanics should be regarded as being 
on an equal footing with standard Hermitian quantum mechanics, and as being completely 
independent of it at the special critical points where the operator Q is singular. 

Having shown that the Pais-Uhlenbeck Hamiltonian defines a physically acceptable 
quantum-mechanical theory, we believe that the techniques and results that we have de- 
scribed here will be of value in quantum field theory. The original motivation of Pais and 
Uhlenbeck was to see if one could avoid the renormalization infinities in theories such as 
quantum electrodynamics by having Feynman propagators that behave as l/k 4 rather than 
as 1/k 2 . The discouraging result that they found was that one could do so but at the price 
of having an energy spectrum without a lower bound. Subsequently, following the devel- 
opment of indefinite-metric theories, it was realized that one could evade this problem and 
have a spectrum that is bounded below, but one would have to pay a different and equally 
unpalatable price, namely, allowing states of negative Dirac norm and violating the physical 
requirement of unitarity. In our work on the Pais-Uhlenbeck model we have been able to 
overcome both the spectral nonpositivity and the norm nonpositivity (ghost) problems. Our 
optimism that ghost problems in quantum field theory can be solved is strengthened by the 
fact that VT techniques have previously been used [U] to show that the quantum-field- 
theoretic Lee model is ghost free. 

There are many possible directions for future research. Perhaps, the most intriguing 
future application of the ideas developed in this paper is to attempt to construct a consistent 
fourth-order derivative quantum theory of gravity in four spacetime dimensions. (Success in 
this endeavor might eliminate the need for ten-dimensional string theory.) Quantum field 
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theory is beyond the scope of this paper, but we note that the Pais-Uhlenbeck model serves 
as a quantum-mechanical prototype for higher- derivative theories, such as conformal gravity 



3l| , which seek to construct a consistent quantum theory of gravity. Roughly speaking, the 



issues involved are illustrated by the fourth-order scalar field theory whose action is 
1 



/ = | / d 4 x 



d^d^d^-iMf + M^d^^ + MfM^ 2 . (101) 



Here, the scalar field can be thought of as representing a typical component of the metric 
fluctuation h^ u around a flat background rj^. The action f llOip gives the equation of motion 

{d 2 - V 2 + M 2 ){d 2 - V 2 + Ml)<P{x, t) = 0. (102) 

With k 2 = k 2 — k 2 , the propagator for this theory is given by 

{d 2 - V 2 + M 2 ){d 2 - V 2 + M 2 )D^{x -y) = 5\x - y), 

D ^ ^ = (k 2 - M 2 )(k 2 - M 2 ) = M 2 - M 2 \k 2 -M 2 ~ k 2 -M 2 ) ' (1 ° 3) 

This propagator exhibits the good 1/k 4 convergence at large k 2 that is needed to make the 
quantum theory renormalizable, but it appears to do so at the expense of having ghost 
states. The identifications u x = (k 2 + Mf) 1 / 2 , u 2 = (k 2 + Mf) 1 / 2 reduce ( TT03D to ([29]), so 
there is hope that there might be a solution to the higher-derivative gravity ghost problem 
that parallels the solution to the Pais-Uhlenbeck-model ghost problem presented here and 
in Ref. [5J. 

To see how things could possibly develop, we note that for the generic higher-derivative 
scalar action I = J d 4 x £(<f), <9 M 0, d^dy/p) the translation invariance of the theory leads to an 
energy-momentum tensor of the form 



T 



dC n ( dC 

771 X 771 



dC 

." + SI ^V> A ~ ^ £ " ( 104 ) 



'9,H \ u( P,v,\, 

This energy-momentum tensor is conserved in field configurations that obey the fourth- 
order-derivative Euler-Lagrange equation 

^-M^) +SA fcr a (105) 

In a phase space formulation of the theory, T^ is replaced by 

T^ u = n^,, + n^^x ~ ri^C. (106) 
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Thus, for the particular action in (110 II) T^ has the form 



2 i\/r2 i2 



Tp, = tv</v + tt^tt\ - - Vflv [7r AK 7r AK - (Mf + M 2 2 )d A 0<9 A + M{ M l 2 



(107) 



and the equation of motion is (I102p . For the metric 77^ = (1,-1,-1,-1), the Hamiltonian 
of the theory is H = J d z x T 00 , where 

Too = vr o + ^vr 2 + 1 -{M\ + M 2 )0 2 - 1 -M\mU 2 - ~tt^ + ^(M 2 + M 2 )0,^\ (108) 



We note that (11081) is similar in structure to -ffpu in fl30l) . and fl 1 8 [) is its covariant 
generalization. Recalling the identifications ixi\ = (k 2 + M 2 ) 1//2 , co> 2 = (k 2 + Mf) 1 / 2 , the 
Hamiltonian (T4"5"l) and commutation relations (Hoi) covariantly generalize to 



# 



d 3 fc 



2(M 2 - M 2 )(P + M 2 )a a a a + 2(M 2 - M 2 )(P + Mi)a 2> - k a 2>k 



+ \l 2 + Ml) 1 / 2 + -{k 2 + Ml) 1 ' 2 



[a a ,a 1]P ] = ^-M^XP + M 2 ) 1 ^]-^-^) 
[a 2>s , a 2jP ] = [2(M 2 - M 2 )(k 2 + M 2 ) 1 / 2 }-^ - k') 
Gh i, Q-9 t' — 0, a-i £, do £/ = 0, Oi e, On p = 0, lai s, a 



n,^ «2,fc'j 



n,fe) "2,/c'J 



n.fe) «2,fc'j 



l,fc; u 2,fc'J 



0. 



(109) 



In (11091) all relative signs are positive. 

Similarly, in the limiting case where M 2 = M 2 = M 2 , the above equations are replaced 
by the covariant generalizations of the Jordan-block Hamiltonian (1541) and the commutation 
relations (15"6T) : 



H 



d s k 



8(F + M 2 f[2b- k b- k + a~ k b- k + b m ] + (P + M 2 ) 1 / 2 ] , 
[a- k , h] = [b- k , ay] = [8(P + M 2 f 2 ]- l 5\k - ¥), 
[oit, %'] = 0, [b k , b k ,} = 0, [a k , b k ,} = 0, [a k , b k ,} = 0. (110) 

Now, there are zero-norm states and this continues to be the case even if we set M 2 = 0. 
Thus, we find zero-norm states in the pure fourth-order conformal gravity case where the 
action can schematically be represented by / = | / d i xd p d v $d il d h '(f) alone 32J. 
C. M. B. is supported by a grant from the U.S. Department of Energy. 
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APPENDIX A: LEHMANN SPECTRAL REPRESENTATION FOR HIGHER- 
DERIVATIVE THEORIES 

To derive the Lehmann spectral representation in a higher-derivative field theory, we 
follow the same procedure as that for a second-order derivative field theory. The derivation 
makes use of three assumptions: the Poincare transformation properties of the interacting 
fields, a completeness relation for the exact energy eigenstates of the interacting theory, and 
the existence of a stable ground state \Q) whose four-momentum we can take to be zero. 

For the case of a self-interacting Hermitian scalar field (j)(x), one introduces the translation 
generator P M = J d 3 xT 0fl , with the scalar field then transforming according to 

[0(x), P„] = id^(x), <f>(x) = e iPx me- lP - x . (Al) 



Given (I A II) . the matrix element of 4>(x) between the vacuum and the one-particle state of 
four- momentum fc" positive k$, and squared mass m^ = k n ■ k n is: 

(n\<j>(x)\k$ = W(0)|A>-^, (k;\<P(x)\n) = (A£W0)|n)e""". (A2) 

Let us provisionally take the completeness relation for the four-momentum eigenstates to 
be of the conventional Dirac form 

$»<n| = l. (A3) 

n 

We can then write the two-point function of the field 4>(x) in the form 

(m^vm = £ i w(o)iA£)iv*"-<--*>. (A4) 

n 

Introducing the spectral function p(q 2 ) defined by 

p(q 2 ) = (27r) 3 E^-^)l^l^°)l^>| 2 %o), (A5) 

n 

we write the two-point function as 

(n|0(aMy)|n) = jH dm 2 p{m 2 ) J ^9(q )5(q 2 - m 2 )^-*). (A6) 

Repeating the same analysis for the (Q\<f)(y)<j)(x)\Q) two-point function, we obtain the 
usual Lehmann representation 

A^(x -y)= i(n\T[<f>(x)<f>{y)]\n) = i(0| [9(x - y )<f>(x)<f>{y) + 6(y - x Q )<f>{y)<t>{x)\ \Q) 

dm 2 p(m 2 )A^ 2) (x-y;m 2 ), (A7) 
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is|. Here, A& 6 ,/- — ' ' 



^,2)1^ - y;wi 



for the time-ordered product of the interacting fields 

[i(Q\T[(f)(x)(j)(y)] |fi)]ftee is the Feynman propagator for a free scalar field with mass m, namely 

the propagator that obeys the free second-order-theory relations 

d 2 _ V 2 + m 2] A free )(x _ y . ^ = ^ _ y)? Aft^* 2 ; m 2 ) = -^^T (A8) 

(Comparing (1A8|) with the fourth-order propagator in (11031) . we note a relative minus sign 
in Fourier space.) 

Equation ( 1A7I) expresses the exact Feynman propagator of the interacting theory as a 
spectral integral over free second-order Feynman propagators with a continuum of mass 
values. No assumption has been made regarding the order of the equation of motion obeyed 
by the interacting field (f>(x). Nevertheless, the free propagator Afp e 2 Jx — y;m 2 ) in (1A7j) 
still satisfies a second-order differential equation. Independent of the structure or order of 
the interacting field equation, the mass-shell condition associated with the eigenstates of 
the exact four-momentum operator P M is still a second-order condition because of Poincare 
invariance. 

The derivation of (1A7j) is generic, and we can apply it to a field theory having a fourth- 
order field equation. To do so we need to show that in the absence of interactions the 
quantity i(Q\T[(p(x)(j)(y)}\Q) can indeed be identified with the free fourth-order propagator 

D<1) «=M^M|(^W-^W) (A9) 

introduced above. In the free fourth-order case associated with the equation of motion 

(a 2 - V 2 + M 2 )(«9 2 - V 2 + M 2 )0(x,t) = 0, (A10) 

we must evaluate the relevant time derivatives of the two-field T-product. 
For the first time derivative we obtain 

^- [(n\T[<i>(x)<f>(y)]\n)] = (n\ [e(x - yo)m<P(y) + %o - x )</>(y)j>(x)] \n) 

+ {n\6(xo-y )[<f>(x),<l>(y)]\n), (All) 

and, with 5(xq — yo)[4>(x), </>(y)] = 0, we obtain 

j?-[(n\T[<j>{x)<i>(y)]\n)] = (n\ \e(x - y )4>(x)<f>(y) + e(y - x )<i>{y)j>(xj\ \n). (A12) 

OXq l j 
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Analogously, for the second time derivative we obtain 

d 2 



8xq 



[(n\T[<j>(x)<p(y)]\n)] = (n\l9(x -yo)(i>(xMy) + e(yo-xo)(J>(y)(J>(x)\\ty 
+ (Q|5(x o -2/o)[0(x),0(y)]|O). 



(A13) 



However, unlike the second-order case, <p(x) and 4>(x) are not canonical conjugates. Rather, 
they are the covariant analogs of the y and y = i[H, y] = —ix variables of the fourth-order 
Pais-Uhlenbeck oscillator theory. With the y and x operators being commuting variables, the 
fourth-order equal-time commutator 5(xq — yo)[<fi(x), <j>{y)\ must thus vanish. Consequently, 
(1A13P reduces to 

d 2 



OXr 



[(fi|T[0(x)0(y)]|fi)] = (Q\ 6(x - y a )4>{x)(j>{y) + %o - x Q )<j>{y)(t>{x) |fi>. (A14) 



For the third time derivative we obtain 
<9 3 



dxl 



[(Q|T[0(x)0(y)]|O)] = (fl| 



0(x o - y ) 



dXn 



c) cj ) (y) + o(yo-xomyf" (r) 



8xq 



+(n\6(x -y )[4>(x),<i>{y)]\n). 



(A15) 



Because <j)(x) and <j)(x) are the covariant analogs of the y and — ix = [H, x] = —ip/'y variables, 
the commutator S(xq — yo)[(j)(x), 4>{y)] also vanishes. Thus, ( 1A15I) reduces to 



a 3 

dXn 



m\T[<j>(x)(t>(ym)} = (Q\ 



d d>(x) 

8{x ~ Vo) o 3 <f>{y) + %o - x Q )(j)(y) 
ox 



Finally, for the fourth time derivative we obtain 
<9 4 



dx A 



[<n|T[^(a;)^(y)]|n>] = (0| 



d d>(x) 
0{xo ~ Vo) a 4 <f>{y) + %o - x Q )(f)(y) 



Oxq 



+ (Q\5(x -y ) 



d 3 (j)(x) 
dx'n 



4{y) 



|n). 



d 3 4>(x) 
dxl 


|n>. 


(A16) 


d A 4>(x) 
8xq 


|n> 




(A17) 



However, d 3 (j)(x)/dx 3 ) is the analog of —ip/'~f = [if,p]/7 = i(a;f+a;|)aT+q , /7, and the commu- 
tator of |/ and g is [?/, g] = i. Thus, the equal-time commutator 5(xq — yo)[d 3 <p(x)/dx 3 ) , <j){y)\ is 
equal to — iS 4 (x — y). Hence, the i times the T-product obeys (11031) . and i(Q\T[(fi(x)(f)(y)]\Q) 
is the fourth-order-theory Green's function we need. 

Note that in the fourth-order case in Fourier space the left side of (1A7I) behaves as 1/fc 4 at 



large k 2 . However, the Fourier transform of A 



free 
(F,2) 






y\ m 2 ) behaves as 1/k 2 . The spectral 
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function p{m 2 ) in (1A5I) is positive definite. Since (1A7j) is a mathematical identity, we thus 
have a contradiction. The difference between the large k 2 behaviors of Afe \s (x — y) and of 
A| r F ee 2 |(x — y) implies that the spectral function p(m 2 ) is not positive definite. Hence, the 
standard Dirac completeness relation in (IA3P for the energy eigenstates cannot be valid in 
the fourth-order case. This analysis immediately generalizes to all higher-order derivative 
theories for which the propagator is even more convergent at large k 2 . We conclude that 
before constructing the Hilbert space appropriate to a higher-derivative theory, we know 
from the outset that the needed inner product cannot be the standard Dirac one. 

A simple modification of the Dirac norm that gives a nonpositive definite p(m 2 ) would 
be to replace (1A3I) by 

5»)(n| = l, (A18) 

n 

where r\ n is ±1, and to replace flA5j) by 

p(q 2 ) = (2nfY,S\K-^\WmK)M%)- (A19) 

n 

However, this choice would violate unitarity. 

When the Hamiltonian is not Hermitian, positivity of the spectral weight function is no 



longer mandatory 33]. However, we need not forego unitarity because we no longer use the 



Dirac inner product. Instead, we use the VT inner product 



me 



m) = 6 nhn , ]T|n>(n|e- Q = l (A20) 



introduced earlier. With such an inner product the spectral function of ( 1A5I) is replaced by 
p(q 2 ) = (2 7 r) 3 ^5 4 (^-g,)(fi|e-^(O)|^)(^|e- s 0(O)|fi)e(g o ). (A21) 

n 

Unlike the spectral function associated with a Dirac norm, this spectral function need 



not be positive definite. To evaluate it, in analogy to (jHj) we set 
d 3 k 

We define right and left vacua according to 



x) = J (27t)3/2 [-zo 1)S e ikl ' x + a 2rk e ik2 ' x - io^e** 1 ** + a^e^] . (A22) 



Oi,fc|0fl) = 0> a 2 - k \0 R ) = 0, (Oi|oi,fc = 0, (0 L |a 2j fc = 0, (A23) 

and in terms of the notation in (IA23P identify 

\Q) = |0 fl >, \k;) = [2(M 2 - M 2 )(k 2 + Mf)^/\- k \0 R ), 

(n\e~ Q = (0 L |, (k;\ = [2(M 2 - M 2 2 ){k 2 + Mf ) 1 /2]V 2 (0 L |a a , (A24) 
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as normalized according to 



(0 L |0 fl > = 1, (*X) = 1, (kl\k 2 u ) = 1. (A25) 



Through use of the commutation relations fll09|) we obtain the on-shell contributions 

wmvixifr-'mM = - (2 ,)3 2( ^ + m?)./ 2( m ? -m|) ' 

(flle-^OJI^XASIe-^OJin) = ^^ + A ^ /1(M? _ M ,y (A») 

to give us precisely the relative minus sign we require. Recalling the minus sign in ( 1A8I) . 
insertion of (1A26|) in the Lehmann representation then yields the propagator (1A9|) . Hence, in 
fourth-order theories compatibility of the Lehmann representation with unitarity is readily 
achievable 
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